Lande g factors and orbital momentum quenching in semiconductor quantum dots 
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We show that the electron and hole Lande g factors in self-assembled III-V quantum dots have a 
rich structure intermediate between that expected for paramagnetic atomic impurities and for bulk 
semiconductors. Strain, dot geometry, and confinement energy significantly modify the effective g 
factors of the semiconductor material from which the dot and barrier are constructed, yet these 
effects are insufficient to explain our results. We find that the quantization of the quantum dot 
electronic states further quenches the orbital angular momentum of the dot states, pushing the 
electron g factor towards 2, even when all the semiconductor constituents of the dot have negative 
g factors. This leads to trends in the dot's electron g factors that are the opposite of those expected 
from the effective g factors of the dot and barrier material. Both electron and hole g factors are 
strongly dependent on the magnetic field orientation; hole g factors for InAs/GaAs quatum dots have 
large positive values along the growth direction and small negative values in-plane. The approximate 
shape of a quantum dot can be determined from measurements of this g factor asymmetry. 
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An individual electron or hole spin in a single semicon- 
ductor quantum dot provides an excellent system for test- 
ing fundamental aspects of quantum dynamics and coher- 
ence. The central quantity characterizing the response of 
an electron or hole spin to an applied magneti^fidd^the 
Lande g factor, has been measured optically 
and electrically 0, H 0, El- For large dots[_ 
such as those defined lithographically!?! or by electrical 
gates 0, the magnetic fields of interest are usually large 
enough that the magnetic length is smaller than the dot 
diameter. In this limit the g factors are closely related to 
those of quantum wells, and theory^ appears to agree 
with experiment 8]. The theoretical situation is much 
less satisfactory for small dots — asymmetric structures 
grown by self assembly in the molecular beam epitaxy 
(MBE) growth process or spherical nanocrystals grown 
by chemical synthesis. Although several phenomena 
known to affect ^-factors in quantum wells llll Il2l have 
been explored in quantum dots[ll El El Til US El, 
the electronic states in dots are discrete, and thus differ 
qualitatively from semiconductors with unbounded mo- 
tion in one or more directions [13. Many dot properties, 
such as the sharply-peaked optical transitions, thus re- 
semble those of atoms more than bulk semiconductors. 
If quantum dots are considered as "artificial atoms" and 
approached with techniques developed for magnetic atom 
dopants in solids |2fJ, then the relevant quantity is the 
ratio of the energy splitting between different angular 
momentum states to the spin-orbit interaction, and for 
a strong confining potential the g factor of dots should 
approach 2. 

Here we show that the g factor of an electron or a 
hole in a quantum dot depends significantly on an atom- 
like property: the quenching of orbital angular momen- 
tum through quantum confinement. To identify this ef- 
fect we must consider also the known bulk-like effects 



on the q fa ctors from dot strain and composition. These 
include [II El El El El El the modification of the elec- 
tron or hole ground-state energy in the dot, the relative 
proportion of dot or barrier material the wavefunction oc- 
cupies, and valence band mixing (typically of heavy and 
light states). In bulk semiconductors the conduction- 
band g factor is|2l|. 
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where Ep is the Kane energy element, E g is the band gap, 
and A is the spin orbit coupling. For unstrained spher- 
ical InAs nanocrystals with hard wall boundary condi- 
tions only E g changes, yet Eq. Q is a very poor predic- 
tor of g factors in these dots (Fig. . For MBE-grown 
InAs/GaAs dots we find that the electron g factor pre- 
dicted from Eq. for the strained InAs is negative, as 
is the known g factor for unstrained GaAs. Yet the elec- 
tron g factors for such dots are positive over almost the 
entire size range (Fig. |2J . Thus the bulk- like approach to 
g factors in these quantum dots, ave ragi ng the g factors 
over the dot and barrier material [l^ LLJ], also fails. The 
competing influence of atom-like and bulk- like effects pre- 
dict that the growth-direction electron g factor increases 
with increasing dot size, whereas considering only bulk- 
like effects leads to the opposite result. Our results agree 
with recent experiments on electron g factors 0. We pre- 
dict hole g factors have large positive values for mag- 
netic fields in the growth direction, and small negative 
values in the in-plane directions(Fig. . and are sensi- 
tive to dot shape. The calculations presented here are for 
InAs spherical nanocrystals and InAs/GaAs MBE-grown 
quantum dots, but our qualitative conclusion — that the 
"artificial atom" viewpoint is vitally important to under- 
standing g factors — applies to all small quantum dots. 
Our results also point the way towards electric-field 
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control of g factors in quantum dots, g factor control via 
manipulation of the electronic wavefunction in quantum 
wells has already been used to control spin precession [22] . 
and to drive spin resonance |23|. Due to the large size 
of quantum dots compared to atoms, moderate voltages 
applied by electrical gates^ can modify the dot shape, 
energy levels, and g, and thus drive spin resonance in a 
static magnetic field. Such control of electron spin reso- 
nance in an individual quantum dot could assist ultrafast 
manipulation of information encoded in electron or hole 
spin, as well as permit single-qubit gate operations for 
quantum conrputationp^) . 

We have calculated quantum dot g factors by calculat- 
ing the spin splittings in a magnetic field of \B\ = 0.1 T. 
The sign of g was determined by examination of the wave 
functions to see if the spin of the lower energy state was 
parallel or anti-parallel to B. The calculations were per- 
formed using 8-band strain dependent k • p theory in 
the envelope approximation with finite differences on a 
real space grid[25j. Material parameters were taken from 
Ref. assuming T = OK. 

The magnetic field was included by coupling to both 
the envelope function and the electron spin. The enve- 
lope was coupled to B by making all difference operators 
covariant using the standard prescription for introducing 
gauge fields on a lattice. For example, 
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where e is the grid spacing and U x (r) is the phase ac- 
quired by an electron hopping from the site at r to the 
site at r + ex. The U's were determined by the re- 
quirement that transport around a plaquette produced 
the Aharonov-Bohm phase corresponding to the encir- 
cled flux, for example 

U x (r)U y {r + ex)U x (r + ey)Ul(f) = cxp(ie 2 B ± e/h)(4) 

where B± is the magnetic field component perpendicular 
to the plaquette. The electron spin was coupled to B 
though a Pauli term for the Bloch functions, given by 
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where /is is the Bohr magneton and a and J are the 
spin matrices for spin 1/2 and 3/2 respectively. The g 



factors for the Bloch functions are 2, |, 



and g for the 



conduction, valence, and spin-orbit bands respectively. 
The Bloch function g factors are determined solely by 
the angular momentum of the Bloch states. 

As these numerical calculations are performed in a fi- 
nite box, numerical artifacts may arise if the total mag- 
netic flux through the box is not an integer number of 



the flux quantum corresponding to a single electron|27|. 
To avoid any such problems, the value of B was modified 
around the edges of the box (within the barrier material) 
so as to make the total flux in each of the x, y, and z 
directions an integer. This allowed the use of both hard 
wall and periodic boundary conditions, and comparison 
between the two boundary conditions was used to estab- 
lish that the box was sufficiently large to avoid finite-size 
artifacts. 

One of the advantages of this approach is that, once 
we settle on a bulk band basis, we do not further trun- 
cate our quantum dot state basis. Because the spin-orbit 
interaction is positive in the bulk semiconductor con- 
stituents of the dot, it is energetically favorable for the 
spin and orbital angular momenta to be anti-parallel to 
each other. (This is the ultimate reason for g < 2). As 
many individual states make positive or negative con- 
tributions to g, a calculation in a truncated dot state 
basis runs the risk of an unbalanced choice of positive or 
negative contributing states. This can produce g > 2, 
unphysical for electrons. 

We now summarize the physical picture we have devel- 
oped to understand the calculations shown in Figs. 11131 
First consider the known origin of the effective g for con- 
duction electrons in a bulk semiconductor [2l|. When a 
magnetic field is applied, the orbital part of the wave- 
function is modified into Landau levels, corresponding to 
quantized orbital angular momentum around the axis of 
the magnetic field. The Zeeman energy now splits the 
lowest Landau level into two spin-polarized Landau lev- 
els, one with spin parallel to the quantized orbital angu- 
lar momentum and one antiparallel. Although the bare 
g = 2 lowers the energy of the parallel spin state and 
raises that of the antiparallel state, the spin-orbit inter- 
action preferentially aligns spin antiparallel to the orbital 
angular momentum. When that effect is absorbed into 
an effective g, it makes g < 2. The situation is modi- 
fied significantly in quantum dots. Instead of beginning 
with a continuous spectrum, which then is modified into 
Landau levels, the spectrum in quantum dots (with or 
without a magnetic field) is discrete. Thus the modi- 
fication of the lowest energy electron state (CI) pair is 
proportional to the ratio of the cyclotron energy to the 
energy splitting between that pair of states and the next 
lowest pair, huj c /(Ec2 — Eci)- For all but the largest 
dots this ratio is very small, and the resulting orbital an- 
gular momentum of the CI states in a magnetic field is 
very small, leading to g — > 2. We refer to this effect as 
quenching of the orbital angular momentum, and note 
that there are similarities between this picture and that 
developed for crystal-field splitting of degenerate d and / 
states for paramagnetic impurities in insulating solids |20j 
and for splitting of the heavy hole and light hole sub- 
band energies in quantum wells|l2j. A significant differ- 
ence, however, is that these previous cases involve the 
splitting of a discrete number of degenerate atomic-like 
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states, whereas the quenching here is of a continuum to 
individual Kramers doublets. 

Figure Q shows our calculations for sp herical InAs 
nanocrystals, and the bulk- like formula|2lJ. The calcu- 
lations were done with a high barrier (10 eV) to avoid 
any leakage of the wavefunction outside the dot. For 
large dots (E g ~ 0.41 eV) the two agree, however they 
rapidly diverge for smaller dots; when the confinement 
energy equals the bulk band gap of InAs the deviation 
from g — 2 predicted by Eq. Q is six times larger than the 
actual value. Hole g factors for these nanocrystals show 
similar evidence of orbital angular momentum quenching 
(not shown). Note that the quenching of orbital angular 
momentum for the electron and hole states is compatible 
with the high fidelity selection rules for generating spin- 
polarized carriers in dots with optical means [2^, because 
the optical transitions connect states with specific angu- 
lar momentum (valence and conduction) whereas the g 
factors probe how much angular momentum admixture 
is possible for B / 0. 

Fig. El shows g factors for the lowest-energy electron 
state and hole state for spherical-cap InAs/GaAs dots 
as the dot size decreases (and E g increases). The dots 
have a height h in the growth direction, [001], and may 
be circular or elliptical (extended in the [110] direction), 
according to the ratio e = d[iio] /dpio] ■ Fig- 121 h and 
e are fixed as the size changes (dots with the same h and 
e and different E g have different size base lengths rfr 110 ] 
and G^To])- For both the electrons and holes the devi- 
ation of g from the quenched case (g = 2) increases as 
Eg increases (and the dot size decreases). For smaller 
dots, excited bound states become squeezed out of the 
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FIG. 1: InAs nanocrystal electron g factor as a 
function of dot size (parameterized by E g ). The 
effects of angular momentum quenching are clearly 
seen as g — > 2 for smaller nanocrystals. 10 eV 
barriers were used to exclude the wavefunction from 
the barriers. The dashed line shows the g fac- 
tor obtained from the bulk formula, using E g for 
the nanocrystal (including the confinement energy). 



dot, and the states with which the ground state mixes 
are increasingly from the continuum of states in the bar- 
rier material. Hence in the limit of a very small dot g 
approaches the value corresponding to the bulk barrier 
material. For very large dots, g should correspond to the 
strained material within the dot. However, this limit re- 
quires an extremely large dot, and is not reached for the 
computationally tractable dot sizes considered here (it is 
in the nanocrystal calculations of Fig. Measurements 
for large dots [a] are shown as well. These measurements 
cannot determine the sign of the g factor; we identify the 
sign to be positive based on our calculations. 

Using Eq.Qwith E g and A replaced by the values cal- 
culated for bulk InAs with the same strain as the InAs 
in the dot leads to the electron g m — 6. Increasing E g to 
include the largest possible confinement energy for a dot 
state gives g w —0.5. The GaAs barrier has g = —0.44, 
so any averaging approach would produce a negative g 
factor, and the g factor should decrease as the quantum 
dot size increases. Due to orbital angular momentum 
quenching the actual electron g factors shown in Fig. [21 
are positive, except for the largest dots considered, which 
have g ~ —0.05, and the g factor increases as the quan- 
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FIG. 2: InAs/GaAs quantum dot g factors for B in the 
[001] direction as a function of the dot size (parameter- 
ized by Eg), for various dot shapes. Solid lines are for 
the lowest energy electron, dashed lines are for the high- 
est energy hole. Experimental values are from Ref. [5j 
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turn dot size increases — the opposite trend. 

As the appproximate dot shape is frequently unknown, 
the [001] hole g factors (Fig. |2J) or in-plane electron g 
factors (Fig. |3J) may be usable as diagnostics for the dot 
height, and the in-plane hole g factors (Fig. [3J for the 
elongation. Use of in-plane electron g factors to diagnose 
dot height may be particularly effective because their de- 
pendence on Eg and e is weak. In-plane hole g factors, 
however, depend strongly on both the height and elon- 
gation, but not on E g . The in-plane hole g factors are 
much smaller in magnitude than the hole g factors along 
the growth direction. This is qualitatively similar to the 
results of Ref. Hfil for Si/Ge quantum dots. However, for 
InAs/GaAs quantum dots the sign of the hole g factor 
differs for the growth direction and in-plane. 

We have shown that the g factors in quantum dots 
are dominated by the atom-like phenomenon of angu- 
lar momentum quenching, leading to electron and hole 
g factor values and trends that differ substantially from 
those seen in bulk semiconductors, and which are much 
closer to the bare value of 2. For InAs/GaAs dots, strain 
and geometric asymmetry also contribute to g, and result 
in strong g factor anisotropy, especially for hole states. 
Careful study of the anisotropy of the hole g factors can 
provide approximate information about the dot shape. 
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FIG. 3: InAs/GaAs quantum dot g factors along the 
[110] and [110] directions as a function of dot size (pa- 
rameterized by Eg). Solid lines are for the lowest energy 
electron, dashed lines are for the highest energy hole. 



Our unexpected results for the electron and hole g fac- 
tors suggest that the "artificial atom" picture promises 
additional surprises for the physics of spin dynamics in 
quantum dots. 

This work was supported by DARPA/ARO DAAD19- 
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